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ABSTRACT 

We have determined beryllium abundances for 25 metal-poor stars based on the high 
resolution and high signal-to- noise ratio spectra from the VLT/UVES database. Our 
results confirm that Be abundances increase with Fc, supporting the global enrich- 
ment of Be in the Galaxy. Oxygen abundances based on [O i] forbidden line implies 
a linear relation with a slope close to one for the Be vs. O trend, which indicates 
that Be is probably produced in a primary process. Some strong evidences are found 
for the intrinsic dispersion of Be abundances at a given metallicity. The deviation of 
HD 132475 and HD 126681 from the general Be vs. Fe and Be vs. O trend favours 
the predictions of the superbubble model, though the possibility that such dispersion 
originates from the inhomogeneous enrichment in Fe and O of the protogalactic gas 
cannot be excluded. 

Key words: stars: abundances - stars: atmospheres - stars: Population II - Galaxy: 
evolution - Galaxy: halo 



1 INTRODUCTION 

The light element beryllium (Be) has very special origins. 
The p rimordial Be abunda nce (on the order of JV(Be/H) = 
lQ- 17 . lThomas et al.|[l99l ) is negligible as predicted by the 
standard Big Bang Nucleosynthesis. Be can not be produced 
by nuclear fusion in the interiors of stars; on the opposite, Be 
will be destroyed by thi s proc ess. It was first proposed by 
I Reeves. Fowler fc Hovlei i|l970l ) that Be can be created by 
spallation reactions between galactic cosmic rays (GCRs) 
and the CNO nuclei in the interstellar medium (ISM). This 
model predicts a quadratic relation between the abundances 
of Be and O (or a slope of 2 in logarithmic plane) , assuming 
that the CNO abundance is proportional to the cumulative 
number of Type II supernovae (SNell) and the cosmic ray 
flux is proportional to the instantaneou s rate of SNe II. How- 
ever, recent observational results (e.g., Gilmore et al.lll~992l ; 



iBoesgaard fc Kindfl993l ; iMolaro et all 



19971 ) showed a lin- 



ear relation between Be and O abundances, which indicates 
that Be may be produced in a prima ry process instead o f 
the standard secondary GCRs process. iDuncan et all (|l997l ) 
suggested that Be can be produced in a reverse spallation of 
C and O nuclei onto protons and a-particles. This process 
will lead to a linear dependence of Be o n O abundances. 
The results from the latest big survey by IBoesgaard et alj 
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(1 19991 . hereafter B99) showed that the slope of the Be vs. 
O trend in logarithmic scale is about 1.5, which makes the 
Be production scenario more complicated and confusing. In 
addition, the exact slope of the Be vs. O relationship de- 
pends on which oxygen indicator is used (see discussion in 
Sect.lHJ. 

If Be is produced in a primary process and the cosmic 
rays were transported globally across the early Galaxy, Be 
abundances should show a very small scatter at a given time. 
This makes the Be abu n dance an ideal c osmic chron ometer 
i|Suzuki fc Yoshiill200lD . iPasauini et ail i|2004 l2007h found 
that Be abundances in globular clusters NGC 6397 and 
NGC 6752 are very similar to that of the field stars with the 
same Fe abundances. Furthermore, the derived ages f rom 
Be abundances based on the model of lValle et al l (2002) arc 
in excellent agreement with the ages determined by main- 
sequence fitting. They suggested that Be is produced in pri- 
mary spallation of cosmic rays acting on a Galactic scale, 
and the refore can be used as a cosmo chronometer. However, 
B99 an d lBoesgaard fc Novickil (|2006^ found strong evidences 
for intrinsic spread of Be abundances at a given metallicity, 
which may indicate that there is al so local enrichment o f 
Be in the Galaxy. But interestingly, IPasauini et al.l (|2005l ) 
found that stars bel onging to the the ac cretion and dissipa- 
tive populations (see lGratton et aT]|2003l for the exact defini- 
tions for these two kinematical classes) are neatly separated 
in the [O/Fe] vs. log(Be/H) diagram, and especially, the ac- 
cretion component shows a large scatter in the [O/Fe] vs. 
log(Be/H) diagram. They proposed that most of the scatter 
in the Be vs. Fe (O) trend may be attributed to the intrinsic 
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spread of Fe and O abundances (probably due to the in- 
homogeneous enrichment in Fe and O of the protogalactic 
gas), rather than Be. 

In this work, we present Be abundances of 25 metal- 
poor stars, for most of which the Be abundances are derived 
for the first time. Oxygen abundances are also determined 
from both [O i] forbidden line and O I triplet lines to inves- 
tigate the chemical evolution of Be with O in the Galaxy. In 
Sect. [2] we briefly describe the observations and data reduc- 
tion. The adopted model atmosphere and stellar parameters 
are discussed in Sect. [3] Sect. [4] deals with the abundances 
determinations and uncertainties. Sect. [5] presents the re- 
sults and discussions, and the conclusions are given in the 
last section. 



2 OBSERVATIONS AND DATA REDUCTION 

Our analysis are based on the high resolution and high 
signal-to-noise ratio spectra of 25 metal-poor main sequence 
and subgiant stars from the archive database of observations 
obtained with UVES, the Ultrav iolet and Visual Echelle 
Spectragraph (|Dekker et alj 120001 ) at the ESO VLT 8 m 
Kueyen telescope. The spectra were obtained during two 
observation runs: April 8-12, 2000 and April 10-12, 2001 
(Programme ID 65.L-0507 and 67.D-0439), both with stan- 
dard Dichroic #1 setting in the blue and red arms. The blue 
arm spectra ranged from 3050 to 3850 A with a resolution 
of 48 000, while the red arm spectra ranged from 4800 to 
6800 A with a resolution of 55 000. 

The spectra were reduced using the standard ESO MIDAS 
package. Reduction procedure includes location of echelle or- 
ders, wavelength calibration, background subtraction, flat- 
field correction, order extraction, and continuum normaliza- 
tion. 



3 MODEL ATMOSPHERES AND STELLAR 
PARAMETERS 

We adopted the one dimensional line-blanketed local 
thermodyn amic equilibrium (LT E) atmospheric model 
MAFAGS jFuhrmann et alj Il99jt ) in our analysis. This 
model utilizes the Kuruczl (fl992'l ODFs but rescales the 
iron abundance by — 0.16dex to matc h the improved so- 
lar ir on abundance of loge(Fe) = 7.51 jAnders fe Grevessel 



Table 1. Stellar parameters adopted in the analysis. 



1989). Individual models for each star were computed with 
a-enhancement of 0.4 dex if [Fc/H] < —0.6 and with the 
mixing length parameter l/H p = 0.5 to get consistent tem- 
peratures from the Balmer lines. 

The effective temperatures were derived by fitting the 
wings of Ha and H/3 lines, and then averaged. iNissen et al.l 
(2002) studied oxygen abundances of a large sample stars, 
which includes all the objects investigated in this work 
(actually our sam ple is a subset of those employed by 
INissen et ail l2002h . They determined the effective tem- 
peratures from the b — y and V — K colour indices 
ba sed on the infrared flux method (IRFM ) calibrations 
of lAlonso. Arribas fe Martinez- Roger] |1996). As shown by 
Fig. [T^, the agreement between the two sets of tempera- 
tures are good for most of the stars wi th a mean dif f erence 
of 15 ± 89 K. For the star G 020-024, INissen et all (|2002T l 



Star 


T e ff 


log 9 


[Fe/H] 




Mass 




K 


cgs 


dex 


km s — 1 


Mq 


liu i uyo^ 


5890 


4.12 


—0 89 


1 2 


0.91 


liu y i OiU 


6030 


4.22 


— 1.20 


1 3 


0.81 


■RF) Q7Q1 fi 
in lj y ( y _lu 


6350 


4.11 


—0.88 


1 5 


1.03 


LLLJ 1UJ ( Zirj 


6005 


4.23 


—0.82 


j 3 


0.87 


± L LJ IV/WU'JO 


5990 


4.43 


— 1.30 


1.2 


0.81 


TTF) 1 1 1 QSO 


5850 


3.94 


— 1 11 


1 2 


0.83 


TTF) 1 1 3(S7Q 
itjj iiju i y 


5740 


3.94 


—0.70 


1 2 


0.94 




5720 


4.40 


—0.73 


j i 


0.80 


un 1 991 Qfl 


5975 


3.85 


— 1.74 


1 5 


0.81 


un 1 9Rfici 


5595 


4.53 


— 1.17 


7 


0.71 


ItJJ !Oi4 1 <J 


5705 


3.79 


— 1.50 


j 4 


0.88 


LLLJ J.'rHJZiOO 


5725 


3.68 


—2.41 


1 5 


0.79 


HD 160617 


5940 


3.80 


-1.78 


1.5 


0.86 


HD 166913 


6050 


4.13 


-1.55 


1.3 


0.77 


HD 175179 


5780 


4.18 


-0.74 


1.0 


0.85 


HD 188510 


5480 


4.42 


-1.67 


0.8 


0.55 


HD 189558 


5670 


3.83 


-1.15 


1.2 


0.95 


HD 195633 


6000 


3.86 


-0.64 


1.4 


1.11 


HD 205650 


5815 


4.52 


-1.13 


1.0 


0.77 


HD 298986 


6085 


4.26 


-1.33 


1.3 


0.81 


CD -30° 18140 


6195 


4.15 


-1.87 


1.5 


0.79 


CD —57° 1633 


5915 


4.23 


-0.91 


1.2 


0.84 


G 013-009 


6270 


3.91 


-2.28 


1.5 


0.80 


G 020-024 


6190 


3.90 


-1.92 


1.5 


0.83 


G 183-011 


6190 


4.09 


-2.08 


1.5 


0.77 



gave a temper ature of 6440 K , whic h is 250 K higher than 
ours. Recently, lAsplund et al.l (|2006T l derived a temperature 
of 6247 K for G 02 0-024 based on Ha p rofile fitting, which is 
very close to ours. INissen et alj (l2002f) likely o verestimated 
the reddening of this star l Asplund et alj 12006). 

The surface gravities were determined from the funda- 
mental relation 
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and 



M bol = V + BC + 5logTT + 5 



(1) 



(2) 



where M is the stellar mass, Mboi is the absolute bolometric 
magnitude, V is the visual magnitude, BC is the bolometric 
correction, and 7r is the parallax. 

The absolute visual mag nitude wer e directly derived 
from the Hipparcos parallax (IESAII1997I ) if available with 
a relative error smaller than 30 per cent. For two stars 
G 020-024 and G 183-011, their uncertainties in parallaxes 
are larger than 40 per cent, so only photometric absolute 
visual magnitude can be adop ted. For G 183-011 we fol- 
lowed the M v , p hot = 4.08 from INissen et al.l (|2002T l. While 
for G 020-024 a large dif ference exists between the re- 
sult of INissen et al.l (|2002T ) (M v , p hot = 4.33) and that of 
lAsplund et ail (|2006h (M v , Dhot = 3. 72). We adopted the 
latter value because lAsplund et al.l ([20061 ) used the spec- 
troscopic Ha index which provides more precise estimate of 
interstellar r eddening excess tha n the photometric H/3 index 
employed by INissen et alJd2002h. The bolometric correction 
was taken from lAlonso. Arribas fe Martine~ Rogerl (|l995l ) 
and the stellar mass was estimated by comparing its posi- 
tion in the log(L/LQ)-logT e a diagram with the evolutionary 
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Figure 1. Comparison of effective temperature, surface gravity and iron abundance with lNissen et al ] (|2002h . 



Table 2. Fe II lines used to determine the iron abundances. 
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tracks o f lYi. Kim fe Demaraud {2003) . The final log g values 
are given in Tabic [T] Our r esults are 0.03 ± 0.11 dex lower 
than those of iNissen et al.l (|2002T l on average. Good agree- 
ment holds on for majority of the stars except for G 020-024, 
whereas a difference of 0.43 dex exists. This is mainly due 
to the very differen t absolute visual magnitude adopted by 
INissen et al.l ((2002) and us as discussed above. The differ- 
ence in Mv (0.61 mag) alone will introduce a difference of 
0.24 dex to the surface gravity according to Equation (1). 

Iron abundances were determined from 12 unblended 
Fen lines with spectra synthesis method. We adopted the 
differ ential logg.f values w ith respect to loge(Fe)Q = 7.51 
from iKorn. Shi fe Gehrenl ll2003h and the van der Waals 
damping constants from lAnstee fe O'Maral l|l99ll . Il995h . 
Our final iron abunda n ces ar e in excellent agreement with 
those of INissen et all (|2002T ), who also derived the iron 
abundances from Fe II lines. The mean difference is only 
—0.02 ± 0.05 dex. The microturbulent velocities were deter- 
mined by requiring that the derived [Fe/H] are independent 
of equivalent widths. 

The typical error for our effective temperature is about 
±80 K. The uncertainty of parallax contributes most to 
the error of the surface gravity. The typical relative er- 
ror of ±15 per cent in parallax corresponds to an error of 
±0.13 dex. In addition, the estimated error of ±0.05 Mq in 
stellar mass translates to an error of ±0.02 dex, while errors 
of ±80 K in effective temperature and ±0.05 mag in BC each 



leads to an uncertainty of ±0.02 dex. So the total error of 
log <j is about ±0.15 dex. It is already noted that the iron 
abundance is insensitive to the effective temperature. The 
uncertainty of [Fe/H] is dominated by the error of surface 
gravity. A typical error of ±0.15 dex in log<? results in an 
error of about ±0.07 dex in [Fe/H]. Combined with the line- 
to-line scatter of about ±0.03 dex, the total error of [Fe/H] 
is about ±0.08 dex. And the error for the microturbulent 
velocity is estimated to be about ±0.2 km s . 



4 ABUNDANCES AND UNCERTAINTIES 
4.1 Oxygen 

As Be is mainly produced by the spallation of CNO nuclei, 
oxygen abundances can be a preferred alternative to iron 
in investigating the galactic evolution of Be. It is also im- 
portant to know the O content when determining Be abun- 
dances, because there are many OH lines around the Bell 
doublet. Therefore, the oxygen abundances were firstly in- 
vestigated here. 

There are several indicators for oxygen abundance: the 
ultraviolet OH, the [O i] 6300 and 6363 A, the infrared 
O I 7774 A triplet and the infrared vibration-rotation OH 
lines. However, OH molecules and O I atoms in the lower 
state of the 7774 A transitions are minority species com- 
pared to total number of oxygen atoms, thus are very sensi- 
tive to the adopted stellar parameters, such as the effec- 
tive temperature and surface gravity. Moreover, line for- 
mations are far from LTE for either the ultraviolet OH 
jfflnkle fe Lambertlll975l ; lAsplund fe Garcia Pere3l200ll ) or 
the O I 7774 A triplet lines (iKiselmarfeOOll - ). In contrast. [O i] 
lines are formed very close to LTE and nearly all the oxy- 
gen atoms in the photosphere of dwarf and giant stars are in 
the ground configuration whi ch provide the low er and upper 
level of the forbidden lines l|Nissen et al.l [2002T ) . Therefore, 
it is believed that the [O i] line is the most reliable indicator 
for oxygen abundances, but the difficulty is that the [O i] 
lines are very weak in dwarf and subgiant stars. 

High resolution and high sinal-to-noise ratio spectra 
covering the infrared O I triplet for 10 stars were avail- 
able from the archived VLT/UVES spectra database. We 
re-reduced the spectra and measured the equivalent widths. 
Since for six of these ten stars O abundances from the O I 
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EW (this work) 

Figure 2. Comparison of Ol triplet equ ivalent widths be- 
tween this work and the literatures {circles: |jmiseh^ i et_alj|2005|; 
diam onds: iNissen fc Schuster! 1 19971: squares: iBoesgaarcTfcKirigl 
Il993t) . 



triplet were previously derived by other authors, we com- 
pare our measurements with those found in literature in 
Fig. It can be seen that the agreement is good on the 
whole. For the rest 15 stars, we collected their Ol 7774 A 
triplet equivalent widths from the literatures directly. Oxy- 
gen abundances were computed with the log gf = 0.369, 
0.223, and 0.002 from lWiese. Fuhr fc Deters! (| 19961 ) in LTE 
first. Then n on-LTE correc tions were applied according to 
the results of iTakedal (I2003T). 

In addition, INissen et all l|2002f ) measured the equiva- 
lent widths of [O i] 6300 A line for 18 main-sequence and 
subgiant stars, of which 15 stars are included in our sam- 
ple. They performed the measurement in a very careful 
fashion, including removing the possible blending telluric 
O2 and H2O lines with the observed spectra of rapidly ro- 
tating B-type stars and subtracting the equivalent width 
of the blending Nil line at 6300.339 A. The typical er- 
ror for the equivalent width of [0 1] 6300 A line was esti- 
mated to be about only ±0.3 mA. For these 15 stars, we 
also determined their oxygen abunda nces with the equiv- 
alent widths of [O 1] 6300 A line from INissen et all <|2002h 
using the accurate oscillator strengt h log o f = —9.72 from 
lAllende Prieto. Lambert fc Asplundl l|200ll ). 

Finally, we derived the oxygen abundances with O I 
7774 A triplet for all the 25 sample stars and with [0 1] 
6300 A line for 15 stars, which are given in Table [3] (the 
reference solar O abundance is loge(O) = 8.77 computed 
from [Ol] 6300 A line using the equivalent width of 4.1 mA 
from INissen et all I2002T) . 

Oxygen abundance from the weak [0 1] 6300 A line is 
not sensitive to stellar parameters. Its uncertainty is domi- 
nated by the error in equivalent width. Normally, a typical 
error of ±0.3 mA in equivalent width corresponds on aver- 
age to an error of ±0.1 dex in oxygen abundance. For the 
infrared Ol triplet, errors of ±80 K in effective tempera- 
ture and ±0.15 dex in gravity each translates to an error of 
±0.05 dex in oxygen abundance. The uncertainties in iron 
abundance and microturbulence nearly have no effect on 



0.301 — ' — ' — ' — ' — ' — ' — ' — ' — ' — 1 — ' — ' — ' — ' — ' — ' — ' — ' — ' — 

o o 

o 

I o " 

U o 

I 0-15- o ° c° ~ 

C o 

I 

i 0.00 ° ° 

o 

-mi l ° ................. . 

-2.5 -2.0 -1.5 -1.0 -0.5 

[Fe/H] 

Figure 3. Comparison of O abundances from [Ol] 6300 A line 
with those from O I 7774 A lines. 

[O/Fe]. A typical error of ±3mA in equivalent width corre- 
sponds to an error of ±0.04 dex. In total, the error of [O/Fe] 
from Ol triplet is around ±0.08 dex. The uncertainty stated 
above is a random one, and the systematic error can be much 
higher, as can be seen from the differences of O abundances 
between [0 1] forbidden line and O I triplet lines discussed 
below. 

As there are 15 stars with oxygen abundances from both 
[O 1] 6300 A and O I 7774 A lines, it is interesting to inves- 
tigate whether these two indicators give consistent oxygen 
abundances. Fig. [3] shows the differences of O abundances 
based on [Ol] 6300 A and Ol 7774 A lines. We find that, 
on average, O abundances from O I 7774 A lines with NLTE 
corrections are 0.14 ± 0.10 dex higher than those from [O 1] 
6300 A line for the 15 sample stars, which mea ns that these 
two in dicators are not consistent in our study. INissen et al.l 
(2002) found that the mean difference between O abundance 
from the forbidden and permitted lines is only 0.03 dex for 
five of their stars, and they concluded that these two in- 
dicators produce consistent O abundances . As a test, we 
reanalyse the O abundances of lNissen et al.l l|2002l )'s sample 
using their stellar parameters and equivalent widths. The 
only differences are the adopted model atmospheres and the 
NLTE corrections for the Ol triplet lines. We found that, 
for the forbid den line, our a bund ances are almost the same 
as those from INissen et all {2002), but for the triplet lines, 
our m ean LTE O abunda nce is about 0.06 dex larger than 
that of lNissen et~al] (120021), and the NLTE correction for the 
triplet lines from lTakedal {2003! ) is 0.06 dex lower than theirs. 
These two factors lead to a total difference of 0.12 dex for the 
permitted lines. Therefore, the different model atmospheres 
(MAFAGS vs. MARCS) and NL TE corrections are respon - 
sible for the differences. Recently. iGarcia Perez et al.l (2006) 
determined O abundances for 13 subgiant stars with [Ol], 
1 and OH l i nes. T hey followed exactly the same method as 
INissen et all l|2002t) . but their results showed that O abun- 
dances based on Ol triplet are on average 0.19 ± 0.22 dex 
higher than that from the forbidden line. 

4.2 Beryllium 

Beryllium abundances were derived by spectra synthe- 
sis of the Ben 3130 A resonance doublet region using 
the idl/fortran siu software package of Reetz (1993). 
It is well known that this spectral region is rich with 
atomic and molecular lines for solar-type stars, which re- 
sults in substantial line absorption and a deficit of con- 
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Figure 4. Spectral synthesis of the Bell doublet region for the 
KPNO Solar Flux Atlas data. 



tinuum. We firstly computed the synthetic solar spectrum 
around Be II doublet region b ased on the l i ne list care- 
fully compiled and tested by IPrimas et alj |l997), and 
then compared them with the integrated solar flux at- 
las of iKurucz et al.l (|l984T ). Some changes were made in 
order to make the theoretical solar spectrum match the 
IKurucz et ail d 1984) solar fl u x atla s best. The major change 
made to the IPrimas et al.l (|l997l ) line list is that we in- 
creased the log 3/ of the Mnll 3131.017 A line by 1.72 dex 
instead of adding a predicted Fel line at 3131.043 A. Similar 
adjust ment was adopted by iKing, Delivannis fc B oesgaard 
(1 19971 ). Based on this adjusted line list, we reproduced 
the solar flux atlas best with A(Be) = 1.12, which is 
in good agreemen t with the result of A(Be) = 1.15 ± 
0.20 derived by [chmielewski, Brault fc Muellerl (|l975h . 
iBalachandran fc Bell ( 1998h found that the standard UV 
continuous opacity of the sun need to be multiplied by a 
factor of 1.6 in order to get consistent oxygen abundances 
from the UV and IR OH lines. With the increased UV con- 
tinuous opacity, they determined the solar Be abundance 
to be 1.40, which is very close t o the meteoritic value 1.42. 
iBell. Balachandran fc Bautistal (120011 ) proposed later that 
the 'missing' UV opacity could be accounted for by the Fe I 
bound-free transitions. However, until now there is no con- 
firmed evidence about this. But one should keep in mind the 
'missing' UV opacity problem. If it does exist, the Be vs. Fe 
(O) trend of this work and all the previous work based on 
the standard UV opacity might change. 

In addition to some strong OH lines, a strong Ti 11 line 
at 3130.810 A also presents in the Bell region. In order 
to minimize its effect on the beryllium abundance as well 
as to provide a constraint on the location of continuum, 
we derived the Ti abundances for our sample stars from 
the Til 5866.461, 6258.110, and 6261.106 A lines. The os- 
cillator strengths for these lines are differentially adjusted 
to produce the solar Ti abundance loge(Ti) = 4.94. For 
five very metal-poor stars with [Fe/H] < —1.8, the Til 
lines are too weak to be detect ed, thus a common value 
of [Ti/Fe] = 0.35 (|Magainl fl989h was adopted. As a mat- 
ter of fact, due to the metal deficiency of these stars, the 
line blending and continuum normalization were much less 
problematic than solar-type stars. The abundances of other 
elements, which are less critical for Be abundance deter- 



mination, were adopted by scaling a solar composition. 
Beryllium abundances were then determined by varying the 
value of Be to best fit the observed line pro fi les. I t was 
reported by Garcia Lopez. Severing fc Gomezl (1995) and 
iKiselman fc Carlssonl l |l996h that non-LTE effects for Bell 
doublet are insignificant, normally smaller than 0.1 dex. We 
took the weaker 3131.066 A line as the primary abundance 
indicator because it's less blended compared to the stronger 
3130.421 A component. 

The uncertainties of Be abundances were estimated 
from the errors of stellar parameters and pseudo-continuum 
location. An error of ±0.15 dex in surface gravity implies 
uncertainties of ±(0.06-0.09) dex, while an uncertainty of 
±0.08 dex in [Fe/H] corresponds to an error of ±0.06 dex. 
Errors due to effective temperature and microturbulence are 
always within ±0.04 dex in total. The error in continuum 
location was estimated to be less than five per cent in the 
worst case, which results in an error of ±0.15 dex at most. 
The final errors for each star are given in Table [3] 

Be abundances for several stars of our sample have been 
published by other authors and they are summarized in Ta- 
ble |31 We can see that, though relatively large scatter ex- 
ists, Be abundances from different researches are not differ- 
ent within uncertainties for majority of the stars. The ex- 
ceptions are HD 160617 and HD 189558. Our Be abundance 
of HP 160617 is 0.5 dex higher than that of iMolaro et ail 
l| 19971 ). This difference is m ostly due to t h e diffe rent stel- 
lar parameters adopted by IMolaro et al.l (Il997l ) and us. 
Their effective temperature and surface gravity are 276 K 
and 0.51 dex lower than ours, respectively, which result in a 
much lower Be abundance. Another exception is HD 189558, 
where a difference of 0.37 dex in Be abund ance exists be- 
tween the result of lBoesgaard fc Kind i| 19931 ) and ours. The 
slight differences in the adopted s tellar parameters c ould 
not produce such a large difference. iRebolo et alj (|l988h de- 
rived a Be abundance of logiV(Be/H) = —12.0 ± 0.4 dex 
for this star with si milar stellar parameters. It is 1 dex 
lowe r than the result of Boesgaard fc King! l|l993l ) . We noted 
that iBoesgaard fc Kind (|l993F determined Be abundances 
by measuring the equivalent width of Bell doublet, which is 
very sensitive to the location of the continuum. It is proba- 
bly that they overestimated t he con tinu um. Moreover, both 
the sp ectra of IMolaro et al.l (|l997l ) and IBoesgaard fc King! 
(1993) were obtained with 3.6 m telescopes. The signal-to- 
noise ratios around Bell region of their spectra are much 
lower than that of this work. 



4.3 Lithium 

It is well known that beryllium can be destroyed in stars by 
fusion reactions at a relatively low temperature (about 3.5 
million K). In order to avoid contamination of our sample 
from the effect of depletion process in stars, we also deter- 
mined the 7 Li abundances for our sample stars. Because the 
destruction temperature for 7 Li is lower than that of Be, if 
r Li is not depleted, Be should not be depleted eithefl 



1 For subgiant stars, Li and Be can be diluted due to the en- 
largement of the convection zone. In this case, Li and Be will 
be diluted by the same percentage. Nevertheless, Be cannot be 
depleted more than Li in any case. 
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Figure 5. Spectrum synthesis for six representative stars. The dots are the observational data, the solid line is the best-fit synthesis, 
and the dotted lines are synthetic spectra with Be abundances of ±0.2 dex relative to the best fit. 



We adopted the oscillator strengths from the NIST data 
base, namely log gf = 0.002 and -0.299 for the Li I 6707.76 
and 6707.91 A lines, respectively. Collisional broadening pa- 
rameters describing van der Waa ls interaction with hydro - 
gen atoms were take n from lBarklem. O'Mara fc Rossi (Il998l ) 
fsee IShi et al.|[2007l in detail). Li abundance s were derived 
by sp ectra synthesis i n LTE . Results from lAsplund et al.l 
l|2006h and IShi et all 120071 ) showed that non-LTE effects 
for Li I resonance lines are insignificant. The typical error 
for our Li abundance was estimated to be about 0.1 dex. 



5 RESULTS AND DISCUSSION 

5.1 Be vs. Fe and Be vs. O relation 

Encouraged by the agreement of Be abundances between 
our results and literatures, we now turn to investigate the 
chemical evolution of beryllium in the Galaxy. As we have 
mentioned before, it is first necessary to investigate whether 
some of our stars are depleted in Be. Fig. [6] displays the 




1.0 - 



-2.5 -2.0 -1.5 -1.0 -0.5 

[Fe/H] 

Figure 6. Li abundances as a function of [Fe/H]. The filled cir- 
cles are the determined Li abundances, while the inverse triangles 
represent the upper limit. 
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Table 3. Abundances and population membership. 



Star 


[Fe/H] 


6300 


m /hi 

777 '4 I,TF-° 

1 1 11 III ±_J 


7774 n T,TF, 


[Ti/Fe] 




A(Be) 


cr(Be) 


Pop b 


HD 76932 


-0 


89 


-0.46 


-0.24( 5 > 


-0.36 





24 


2.00 


0.73 





1 1 





HD 97320 


— 1 


20 


-0.85 


-0.56( 6 > 


-0.66 





21 


2.32 


0.43 





17 





HD 97916 


— 


88 




-0.06W 


-0.27 





23 


< 1.23 


< -0.76 






1 


HD 103723 


-0 


82 


-0.68 


-0.33( 6 > 


-0.44 





15 


2.22 


0.51 





17 


1 


HD 106038 


— 1 


30 




-0.62( 3 ) 


-0.70 





21 


2.55 


1.37 





12 


1 


HD 111980 


— 1 


11 


-0.76 


-0.38( 6 > 


-0.50 





29 


2.31 


0.67 





13 


1 


HD 113679 


— 


70 


-0.43 


-0.15( 6 > 


-0.28 





32 


2.05 


0.87 





12 


1 


HD 121004 


— 


73 


-0.42 


-0.19( 2 ) 


-0.27 





28 


< 1.18 


0.94 





15 


1 


HD 122196 


— 1 


71 




-1.07( 6 ) 


-1.17 





28 


2.28 


-0.51 





1 1 





HD 126681 


— 1 


17 


-0.72 


-0.65( 4 > 


-0.70 





30 


1.48 


0.90 





12 





HD 132475 


-1 


50 


-1.09 


-0.82( 4 ) 


-0.91 





27 


2.23 


0.62 





13 


1 


HD 140283 


-2 


il 


-1.61 


-1.65( 3 ) 


-1.73 






2.16 


-0.94 





14 


1 


HD 160617 


-1 


78 


-1.34 


-1.24( 3 ) 


-1.33 





23 


2.25 


-0.41 





12 


1 


HD 1 6691 3 


-1 


55 


— 1.16 


— 0.8l( 4 ' 


—0.90 





30 


2.32 


0.27 





14 


o 


HD 175179 


-0 


74 




-0.13( 6 ) 


-0.24 





32 


< 0.87 


0.73 





15 





HD 188510 


-1 


67 




-0.98( 5 ) 


-1.02 





31 


1.48 


-0.25 





13 





HD 189558 


-1 


15 


-0.73 


-0.54( 1 > 


-0.64 





25 


2.24 


0.64 





14 





HD 195633 


-0 


64 


-0.55 


-0.16( 6 ) 


-0.34 





06 


2.25 


0.53 





18 


2 


HD 205650 


-1 


13 


-0.69 


-0.48( 4 > 


-0.54 





21 


1.70 


0.51 





19 





HD 298986 


-1 


33 


-0.93 


-0.70( 6 ) 


-0.79 





15 


2.26 


-0.04 





12 


1 


CD-30°18140 


-1 


87 




-1.09(3) 


-1.18 






2.21 


-0.35 





15 


1 


CD-57°1633 


-0 


91 




-0.5l( 6 ) 


-0.60 





01 


2.15 


0.31 





18 


1 


G 013-009 


-2 


28 




-1.54(3) 


-1.65 






2.21 


-0.84 





13 


1 


G 020-024 


-1 


92 




-1.19(3) 


-1.30 






2.19 


-0.72 





17 


1 


G 183-011 


-2 


08 




-1.27< 6 ) 


-1.36 






2.21 


-0.61 





11 


1 



a Sou rces o f equivalent width for O I t riplet: (l)ICavallo et al.l dl997f) . (2l lNissen fc Schuster] (I1997T ), (31 iNissen et al.l 
||2002| 1, (4l lGratton et alj d2003h ■ (51 Ijonsell et al.l d2005l l. (61 measured from archival UVES spectra (Programme 
ID 65.L-0507). 

b Population membership: - dissipative component; 1 - accretion component; 2 - thin disc. 



Table 4. Comparison of Be abundances with those from literatures. 



Star (1) (2) (3) (4) (5) This work 



HD 76932 


-11.04 ±0.11 


-11.45 ±0.18 


-11.21 ±0.21 


-11.17±0.05 




-11.27 ±0.14 


HD 132475 










-11.43 ± 0.12 


-11.38 ±0.13 


HD 140283 


-12.78 ±0.14 


-13.07 ±0.20 


-12.91 ±0.17 


-13.08 ±0.09 




-12.94 ±0.14 


HD 160617 






-12.90 ±0.27 






-12.41 ±0.12 


HD 166913 






-11.77 ±0.15 






-11.73 ±0.14 


HD 189558 


-10.99 ±0.15 










-11.36 ±0.14 


HD 195633 


-11.21 ±0.07 








-11.34 ±0.11 


-11.47 ±0.18 



References: (11 IBoesgaard fc King! l ll993ll . (2) iGarci'a Lopez et alj l ll995ll . (3) iMolaro et all l ll997ll , (4) B99, (5) 
iBoesgaard fc Novickil ||2006|1 . 



Li abundances as a function of metallicity. We can see that 
six stars in our sample are obviously depleted in 7 Li and 
another two seem to be slightly depleted, while the star 
HD 106038 has an exceptionally high 7 Li abundance, about 
0.3 dex higher than the Spite plateau (it also has an ab- 
normally high Be abundance, see discussion in Sect. 15.21 for 
this star). Among those stars with depleted 7 Li, only one 
(HD 97916, denoted by filled inverse triangle in Fig. [7] and 
Fig. [8} is also depleted seriously in Be, while the rest (de- 
noted by filled squares in Fig. [7] and Fig. [5} seem to have 
normal Be abundances at their metallicities. 

Excluding HD 106038 and HD 97916 (abnormally high 
Be abundance and seriously depleted in Be, respectively; 
these two stars will not be included in the analysis of Be vs. 



O trend either), the relation between Be and Fe abundances 
can be well represented by linear fitting 

log AT(Be/H) = (1.15 ± 0.07) [Fe/H] - (10.24 ± 0.10) 

One may note that HD 132475 and HD 126681 seem to de- 
viate from the trend far beyond the uncertainties. If we ex- 
clude these two stars (see discussion in Sect. [5T2]| . the rela- 
tionship will be 

log N(Be/R) = (1.10 ± 0.07) [Fe/H] - (10.37 ± 0.10) 

Our result is in reasonable agreement with the result 
log AT(Be/H) = (0.96 ± 0.04) [Fe/H] - (10.59 ± 0.03) of B99, 
considering the relatively smaller metallicity range of our 
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Figure 7. Be abundances against [Fc/H] [filled circles: stars 
without 7 Li depletion; filled squares: stars with depleted 7 Li; filled 
inverse triangle: upper limit Be abundance for HD 97916; filled 
star: solar meteoritic Be abundance; open circles: data from B99). 
The solid line is the best linear fitting (not including HD 106038, 
HD 97916, HD 132475 and HD 126681) with a slope of 1.1 for our 
Be vs. Fe trend, while the dotted lines represents the best fitting 
for the B99 data. 



sample. The overall increase of Be with Fe suggests that Be 
is enriched globally in the Galaxy. 

As the yields of Be is believed to be correlated with 
CNO nuclei directly, it's more meaningful to investigate the 
relationship between Be and O abundances. Due to the in- 
consistent oxygen abundances based on the forbidden and 
triplet lines in our study, it is necessary to investigate their 
relations with Be abundance separately. Fig. [8] shows the 
trend of Be with O abundances. Again, the Be abundances 
increase linearly with increasing [O/H] both for the forbid- 
den line (though with relatively large scatter partly due to 
the small sample number) and triplet lines. The relation- 
ships are best represented by 



log ./V(Be/H) 
logiV(Be/H) 



(1.55 ±0.17) [O/H] 
(1.36 ± 0.09) [O/H] 



(10.29 ±0.15) [Oi] 
(10.69 ±0.08) Oi 



If we exclude HD 132475 and HD 126681, the relationships 
will be 



log AT(Be/H) 
logiV(Be/H) 



(1.49 ±0.16) [O/H] 
(1.30 ±0.08) [O/H] 



(10.42 ±0.15) 
(10.81 ±0.08) 



[Oi] 
Oi 



Our result based on O I triplet lines is slightly flatter 
than the result log 7Y(Be/H) = (1.45±0.04) [O/H] - (10.69± 
0.04) of B99. This can be partly due to our smaller [Fe/H] 
range as mentioned above. Besides, the O abundances of 
B99 were averaged from the UV OH and infrared O I triplet 
lines. They thought that such a result (a slope of roughly 
1.5 for Be vs. O) is neither consistent with the secondary 
process, nor the primary process. However, the secondary 
process added with some chemical evolution effects, such as 
an outflow of mass from the halo, indicates that there would 



be a quadratic relation only at the very lowest metallicities 
and a progressive shallowing of the slope to disc metallicities, 
for example a slope of 1.5 between [Fe/H] = —2 and — 1 
(B99). So they suggested that this process is most consistent 
with their results. 

As discussed above, [O i] forbidden line is the most 
reliable O abundance indicator for its insensitivity to the 
adopted stellar parameters as well as the non-LTE effects. 
Our oxygen abundances from [O i] 6300 A line produces a 
'moderate' s l ope 1 .49 for the Be vs. O trend. However, 
iNissen et alj l|2002h studied the effects of 3D model atmo- 
spheres on the derived O abundance from [O i] forbidden 
lines. They found that O abundances based on [O i] will de- 
crease if 3D models are applied. Especially, the decreasing 
amplitude inc reases with decrea sing metallicity (see fig. 6 
and table 6 of lNissen et~alll2002h . W hile 3D effects o n Be n 
doublet are negligible according to iPrimas et all (2000). 
This means that the slope will be closer to one than our 
present result when the 3D effects were considered for the 
[Oi] forbidden lines. This implies that the Be production 
scenario is probably a primary process. 



5.2 Special stars: hints on the Be production 
scenario 

HD 106038 is very special for its exceptionally overabun- 
dance in Li and Be as mentioned before. Its Li abundance 
is 0.3 dex higher than the Li plateau, while its Be abun- 
dance is about 1.2 dex higher than that of normal stars 
with the same metallicity. The Be abundance of this star 
is even similar to the solar meteoritic val u e. Su ch a Be- 
nch star is extrem ely rare. lAsplund et all (|2006T ) derived 
A(Li) = 2.49 and ISmilianic et al.1 i|2008) reported a Be 
abundance of logiV(Be/H) = —10.60 for this star, both in 
good agreement w i th our results. In addition to Li and Be, 
INissen fe Schust er (1997) showed that this star also has ob- 
viously enhanced abundances o f Si, Ni, Y, and Ba. B ased 
on its special abundance pattern, ISmili anic et af] (120081 ) sug- 
gested that HD 106038 is most probably formed in the vicin- 
ity of a hypernova (HNe). 

In addition to HD 106038, another two stars, namely 
HD 132475 and HD 126681, seem to stand out of the Be 
vs. Fe and Be vs. O trends distinctly. Their Be abun- 



dances are about 0.6 and 0.5 dex higher, respectively, 
than that of most stars with the same metallicities. 
iBoesgaard fe Novickil l|2006l ) also found an abnormally high 
Be abundance for HD 132475 (0.5 dex above the normal 
stars at its metallicity). Their sample included another star 
BD±23°3912, which has very similar atmospheric param- 
eters as HD 132475 but very different Be abundance. In 
fact, BD+23°3912 has a Be abundance matching perfectly 
the linear Be vs. Fe relation. Combined with another star 
HD 94028 also with excess Be abundance found by B99, 
IBoesgaard fe Novickil (|2006h concluded that dispersion in Be 
abundances does exist at a given metallicity, which implies 
a local enrich ment of Be in the Gal axy. 

However. iPasquini et al.l (|2005l ) proposed that such dis- 
persion could be mostly attributed to the scatter of Fe 
and O abundances, rather than Be, as we ha ve men- 
tioned in Sect. \T\ Following IPasquini et all (|2005l ), we cal- 
cul ated the space velocit i es usi ng the method presented 
by IJohnson fe Soderbloml (| 19871 ). and determined the or- 
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[O/H] [O/H] 

Figure 8. Be abundances as a function of [O/H], a) Results from [Ol] forbidden lines, b) Results from Ol permitted lines. The symbols 
and lines have the same meaning as Fig. [7] Remember that O abundances of B99, for both of the two panels, are the mean abundances 
based on UV OH and infrared Ol triplet lines. 



bital parameters based on the Galactic mass model of 
I Allen fc Santillan|[|l99ll ) for our sample stars. Input parame- 
ters, such as radial velocities, parallaxes and proper motions 
were obtai ned from the SIMBAD database. According to the 
criteria of lGratton et al.l (|2003l ). fifteen stars in our sample 
belong to the accretion component, nine stars belong to the 
dissipative component and one star belongs to the thin disc. 

Fig. H] shows [O/Fe] vs. [Fe/H] and [O/Fe] vs. 
log A f (Be/H) for our sample stars, based on the results from 
both [O i] forbidden line and O I triplet lines. One may find 
that no clear separation exists between the two components 
in the [O/Fe] vs. [Fe/H] diagram, though the accretion com- 
ponent shows a relatively larger scatter than the dissipative 
component. However, in the [O/Fe] vs. logA r (Be/H) dia- 
gram, the two populations are distinctly different, and espe- 
cially, the accretion component shows a much larger scatter 
compared to the dissipat ive component. Our r esults agree 
well with the findings of IPasauini et all l|2005h . They pro- 
posed that such results support the idea that the formation 
of the two components took place under significantly dif- 
ferent conditions: an inhomogeneous, rapidly evolving 'halo 
phase' for the accretion component, and a more chemically 
homogeneous, slowly evolving 'thick disc phase' for the dis- 
sipative component. The large scatter of the accretion com- 
ponent in the [O/Fe] vs. logA r (Be/H) diagram may reflect 
the inhomogeneous enrichment in oxygen and iron of the 
halo gas. We note that, for our Be-rich stars, HD 106038 
and HD 132475 belong to the accretion component, while 
HD 126681 belongs to the dissipative component. Another 
Be-rich st ar HD 94028, first discovered by B99 and later con- 
firmed by iBoesgaard fc Novickil ([2006), is also classified as 
a dissipative c o mpon ent star according to the definitions of 
iGratton et al.l (|2003l ). While the deviation of the accretion 
component stars could be due to the inhomogeneous enrich- 
ment in Fe and O of the halo gas, HD 126681 and HD 94028 



cannot be interpreted in this way. However, one should keep 
in mind that stellar kinematics is only of statistical meaning 
in describing the Galactic populations. It has been shown by 
many previous work that substantial overlap exists between 
the halo and thick disc stars. Therefore, it is dangerous to 
attribute an individual star to one stellar population and 
accordingly derive a firm conclusion. Moreover, we didn't 
find any distinct differences in the abundance pattern be- 
tween HD126681/HD 94028 and typical halo stars. So the 
possibility that dispersion in Be vs. Fe and Be vs. O trend 
originates from the inhomogeneous enrichment in Fe and O 
of the protogalactic gas cannot be excluded. 

As an alternative, the scatter in Be can be in- 
terpreted by the so-called s u perbu b ble (SB) model 
jHigdon. L ingcnfcltcr & Ramaty 1998]; iParizot fc Drurvl 
1 19991 ; iParizotJ l2000bl ; iRamatv et all l2000l ). The basic idea 
of the SB model is that repeated SNe occurring in an 
OB association can generate a superbubble, in which the 
CNO nuclei (ejected by SNe) mixed with some ambient, 
metal-poor material are accelerated onto the metal deficient 
material in the supershell and at the surface of the adjacent 
molecular cloud, and then broken into smaller atoms like Li, 
Be and B. The produced light elements are then mixed with 
other SNe ejecta as well as the ambient, metal-poor gas. 
However, such a mixing cannot be perfect, and new stars 
can form before all the massive stars explode or the induced 
light elements production occurs. As a result, scatter in 
the abundances o f light elements for a gi ven S B may occur 
jParizotl l2000al ). IBoesgaard fc Novickil (|2006l ) noted that 
Na, Mg, Si, Y, Zr, and Ba abundances of HD 132475 are 
typically 0.2 dex above the mean values of the o ther stars 
at that me tallicity according to the results of iFulbrightl 
l|2000l . l2002l ). We also noted that Y and Ba abundances of 
HD 126681 are roughly 0.2 and 0.15 dex higher, respectively, 
than the mean values of the remaining sample as found by 
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Figure 9. a) [O/Fc] vs. [Fe/H] based on Ol triplet lines, b) [O/Fc] vs. [Fe/H] based on [Ol] forbidden line, c) [O/Fe] vs. log7V(Be/H) 
based on O I triplet lines, d) [O/Fe] vs. log N (Be/H) based on [O i] forbidden line. Open squares represent the accretion component, filled 
squares represent the dissipative component, and the asterisk represents the only thin disc star HD 195633 in our sample. The open 
square with a left pointinq arrow represents the upper limit Be abundance for HD 97916. 



iNissen fe Schusterl |l997h . and the a-elements of this star 
are in the upper range of their sample. Na and a-elements 
are typical ejecta of SNell, and Y, Zr, and Ba, though 
not very efficiently, can also be produced by the r-process 
in SNell. It is probably that stars like HD 132475 and 
HD 126681 were formed from the material that underwent 
enrichment of light elements and SNell ejecta but not 
plenty of dilution process in SBs. In addition, the SB model 
also predicts a primary process for the Be production, 
which is consistent with our Be vs. O trend. 



Except for the possibilities stated above, some other 
scenarios for the overabundance of Be can be excluded. It 
is unlikely that the Be-rich sta rs were accreted from the 
satelli te sy stems of our Galaxy . IShetrone. Cote fe Sargent! 
l|200ll ) and lTolstov et al.l (|2003h found that abundance pat- 
terns among the dwarf spheroidal (dSph) galaxies stars are 
remarkably uniform. We note that a-elements and Y abun- 
dances of the dSph stars are obviously lower than the Be-rich 
stars. The possibility that the overabundance of Be in our 
Be-rich stars could be due to the accretion of a planet or 
planetesimals debris can also be excluded. If the excess Be 
in our Be-rich stars were accreted from some material hav- 
ing similar composition as chondrites meteorites, then the 
mass of the accreted iron would be even larger than the to- 
tal mass of iron in the surface convective zone of the star, 
which is certainly impossible. 



6 CONCLUSIONS 

We have derived Be abundances for 25 main sequence and 
subgiant stars spanning the range —2.5 < [Fe/H] < —0.5. 
Relations between Be and Fe as well as Be and O are investi- 
gated. The Be vs. Fe trend can be well represented by a lin- 
ear relation with a slope of 1.1. This result is in good agree- 
ment with that of B99, and suggests that Be is enriched glob - 
ally in the Galaxy, as proposed bv lSuzuki fe Yoshiil (2001). 
Our Be abundances increase linearly with increasing [O/H] 
based on both the [O i] 6300 A and O I triplet lines, but with 
slightly different slopes. O abundances based on O I triplet 
gives a slope of 1.30 between [Be/H] and [O/H]. This is a lit- 
tle flatter than the result of B99, which may be partly due to 
different metallicity range. The most reliable O abundance 
indicator, [O i] forbidden line gives a steeper relationship (a 
slope of 1.49). However, this slope will decreas e if 3D effects 
are to ok into account according to the results of lNissen et al.l 
(2002), which means that the production process of Be is 
probably a primary process. 

Moreover, we find some strong evidences for the in- 
trinsic dispersion of Be abundances at a given metallicity. 
The special abundance pattern of HD 106038, especially its 
exceptionally high Be abundance, can be interpreted most 
consistently only if the material which formed HD 106038 
was contaminated by t he nucleosynthetic products of a HNe 
l|Smilianic et alj l2008h . The deviations of HD 132475 and 
HD 126681 from the general Be vs. Fe and Be vs. O trend 
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can be interpreted by the SB model. However, the possibil- 
ity that such dispersion originates from the inhomogeneous 
enrichment in Fe and O of the protogalactic gas cannot be 
excluded. 
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